Phase transition and hysteresis in scale-free network traffic 
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We model information traffic on scale-free networks by introducing the node queue length L 
proportional to the node degree and its delivering ability C proportional to L. The simulation gives 
the overall capacity of the traffic system which is quantified by a phase transition from free flow to 
congestion. It is found that the maximal capacity of the system results from the case of the local 
routing coefficient slightly larger than zero, and we provide an analysis for the optimal value of (f). 
In addition, we report for the first time the fundamental diagram of flow against density, in which 
hysteresis is found, and thus we can classify the traffic flow with four states: free flow, saturated 
flow, bistable and jammed. 
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Complex networks can describe many natural and so- 
cial systems in which lots of entities or people are con- 
nected by physical links or some abstract relations. Since 
the discovery of small-world phenomenon by Watts and 
Strogatz appeared in Nature in 1998, and scale-free 
property by Barabasi and Albert 0] one year later in Sci- 
ence, complex networks have attracted growing interest 
among physics community [3j, A, S Hi 0] ■ pointed out 
by Newman, the ultimate goal of studying complex net- 
works is to understand how the network effects influence 
many kinds of dynamical processes taking place upon 
networks [1]. One of the dynamical processes, traffic of 
information or data packets is of great importance to be 
studied for the modern society. Nowadays we rely greatly 
on networks such as communication, transportation, the 
Internet and power systems, and thus ensuring free traf- 
fic flow on these networks is of great significance and re- 
search interest. In the pass several decades, a great num- 
ber of works on the traffic dynamics have been carried out 
for regular and random networks. Since the increasing 
importance of large communication networks with scale- 
free property such as the Internet the traffic flow on 
scale-free networks has drawn more and more attention 
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Researchers have proposed some models to mimic the 
traffic on complex networks by introducing the random 
generation and the routing of packets [1, m M, M M 
[Tsj . Arenas et al. suggest a theoretical measure to 
investigate the phase transition by defining a quantity 
[lOj , so that the state of traffic flow can be classified to the 
free flow state and the jammed state, where the free flow 
state corresponds to the number of created and delivered 
packet are balanced, and the jammed state corresponds 
to the packets accumulate on the network. 
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Many recent studies have focused on two aspects to 
control the congestion and improve the efficiency of trans- 
portation: modifying underlying network structures or 
developingbetter route searching strategies in a large 
network J24[. Due to the high cost of changing the in- 
frastructure, the latter is comparatively preferable. In 
this light, Echenique et al., Wang et al. and Yin et al. 
suggest traffic models based on the local information or 
the local integration of static and dynamic information 
[TtI . [TsI . I13 |. Yan et al. propose a efficient routing 
strategy based on the knowledge of the whole topology 
[20l | . They find that the efficient path results in the redis- 
tributing traffic loads from central nodes to other noncen- 
tral nodes, and the network capability in handling traffic 
flow is improved more than 10 times by optimizing the 
efficient path. 

However, previous studies usually assumed that the 
capacity of each node, i.e., the maximum queue length of 
each node for holding packets, is unlimited and the node 
handling capability, that is the number of data packets 
a node can forward to other nodes each time step, is 
either a constant or proportional to the degree of each 
node. But, obviously, the capacity and delivering ability 
of a node are limited and variates from node to node in 
real systems, and in most cases, these restrictions could 
be very important in triggering congestion in the traffic 
system. 

Since the analysis on the effects of the node capacity 
and delivering ability restrictions on traffic efficiency is 
still missing, we propose a new model for the traffic dy- 
namics of such networks by taking into account the max- 
imum queue length L and handling capacity C of each 
node. The phase transition from free flow to congestion 
is well captured and, for the first time, we introduce the 
fundamental diagram (flux against density) to character- 
ize the overall capacity and efficiency of the networked 
system. Hysteresis in such network traffic is also pro- 
duced. 

To generate the traffic network, our simulation starts 
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with the most general Barabasi- Albert scale-free network 
model which is in good accordance with real observation 
of communication networks In this model, starting 
from mo fully connected nodes, one node with m links 
is added at each time step in such a way that the prob- 
ability Ili of being connected to the existing node i is 
proportional to the degree kt of the node, i.e. 11^ 
where j runs over all existing nodes. 

The capacity of each node is restricted by two param- 
eters: (1) its maximum packet queue length L, which is 
proportional to its degree k (a hub node ordinarily has 
more memory): L ^ a x k; (2) the maximum number 
of packets it can handle per time step: C — /3 x L. 
Here < /? < 1 simply shows that the maximum 
number of handled packets is less than the maximum 
packet queue length L. Motivated by the previous mod- 
els [i,[l0l[lllil,lii,[I3,[ll, the system evolves in parallel 
according to the following rules: 

1. Add Packets - Packets are added with a given rate R 
(packets per time step) at randomly selected nodes and 
each packet is given a random destination. 

2. Navigate Packets - Each node performs a local 
search among its neighbors. If a packet's destination is 
found in its nearest neighborhood, its direction will be 
directly set to the target. Otherwise, its direction will be 
set to a neighboring node h with preferential probabil- 

ity: Ph = ^ . Here the sum runs over the neighboring 

nodes, and is an adjustable routing parameter in that 
the packets are more likely to be forwarded to high de- 
gree nodes when </> > 0. It is assumed that the nodes are 
unaware of the entire network topology and only know 
the neighboring nodes' degree ki. 

3. Deliver Packets - At each step, all nodes can deliver 
at most C packets towards its destinations and FIFO 
(first-in-first-out) queuing discipline is applied at each 
node. When the queue at a selected node is full, the 
node won't accept any more packets and the packet will 
stay at the site and wait for the next opportunity to be 
forwarded. Once a packet arrives at its destination, it 
will be removed from the system. As in other models, we 
treat all nodes as both hosts and routers for generating 
and delivering packets. 

We first simulate the traffic on a network of = 1000 
nodes with mO — m — 5. To characterize the sys- 
tem's overall capacity, we first investigate the incre- 
ment rate rj of the number of packets in the system: 
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ri{R) — limf 



. Here ANp = Np{t + At) - Np{t) 



At ■ ^^^^"^ . 
with (...) takes average over time windows of width At. 

Obviously, ri{R) = corresponds to the cases of free flow 

state, which is attributed to the balance between the 

number of added and removed packets at the same time. 

As R increases, there is a critical Rc at which Np runs 

quickly towards the system's maximum packet number 

and r]{R) increases suddenly from zero, which indicates 

that packets accumulate in the system and congestion 

emerges and diffuses to everywhere. 

Hence, the system's overall capacity can be measured 
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FIG. 1: (color online). The overall capacity of a network with 
iV = 1000, mO = m = 5, Q = 1(a), a = 2(b) and f3 = 0.2. 
The capacity is characterized by the critical value of Rc for 
different <fi. In (a), a = 1, <j)optimai = 0.3 and 7?^""" = 18.7. 
In (b), a = 2, <l)opUm.ai = 0.1 and K^"'' = 42.2. In both 
cases, the maximum of Rc corresponds to a. <j) slightly greater 
than zero marked by a dash line. The data are obtained by 
averaging Rc over 10 network realizations. 



by the critical value of Rc below which the system can 
maintain its efficient functioning. Fig[T] depicts the vari- 
ation of Rc versus (j). The maximum overall capacity 
occurs at (j) slightly greater than 0.0 with i?™"^ = 18.7 
at = 0.3 for a = 1 (a) and = 42.2 at (/) = 0.1 for 

a — 2 (b). The results are averaged from 10 simulations. 

The analytical estimation of Rc is too complicated for 
our routing model. In a recent paper p3| , Germano and 
de Moura present analytical work on the rather simple 
traffic of particle hopping in complex networks. In the 
following, we provide an analysis for the optimal value 
of (p corresponding to the peak value of Rc. In the case 
oi (j) — Q, packets perform random-like walks if the max- 
imum queue length restriction of each node is neglected. 
The random walk process in graph theory has been exten- 
sively studied. A well-known result valid for our analysis 
is that the time the particle spends at a given node is 
proportional to the degree of such node in the limit of 
long times [255. Similarly, in the process of packet de- 
livery, the number of received packets (load) of a given 
node averaging over a period of time is proportional to 
the degree of that node. Note that the packets deliver- 
ing ability C of each node assumed to be proportional 
to its degree, so that the load and delivering ability of 
each node are balanced, which leads to a fact that no 
congestion occurs earlier on some nodes with particular 
degree than on others. Since in our traffic model, an oc- 
currence of congestion at any node will diffuse to the en- 
tire network ultimately, no more easily congested nodes 
brings the maximum network capacity. However, tak- 
ing the maximum queue length restriction into account, 
short queue length of small degree nodes make them more 
easily jammed, so that routing packets preferentially to- 
wards large degree nodes slightly, i.e., <f> slightly larger 
than zero, can induce the maximum capacity of the sys- 
tem. 

This also explain the difference in the position of i?™"^ 
of Fig[l](a) and FiglTJb). Comparing with the case of 
a = 2, the small degree nodes are more easy to jam 
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FIG. 2: (color online). Average travel time for a network with 
N = 1000, mo = m = 5, Q = 2 and /3 = 0.2. (a) Average 
travel time (T) versus (f) for i? = 1, 2 and 5. The data are 
truncated because the system jams when is either too large 
or too small, (b) The variation of (T) versus R when <p is 
fixed. The data are also truncated when the system jams. 



when a = 1, so a greater </> is needed to achieve a more 
efficient functioning of the system. One can also conclude 
that the optimal will be zero if a is large enough. 

Then we simulate the packets' travel time which is 
also important for measuring the system's efficiency. In 
FiglUa), we show the average travel time (T) versus 
under the conditions of i? = 1, 2 and 5. In the free- flow 
state, almost no congestion on nodes occurs and the time 
for packets waiting in the queue is negligible, therefore, 
the packets' travel time is approximately equal to their 
actual path length in map. But when the system ap- 
proaches a jammed state, the travel time will increase 
rapidly. One can see that when (j) is slightly greater than 
zero, the minimum travel time is obtained. In Fig[21[b), 
the average travel time is much longer when is nega- 
tive than it is positive. These results are consistent with 
the above analysis that a maximum Rc occurs when is 
slightly greater than zero. Or, in other words, this can 
also be explained as: when (j) > 0, packets are more likely 
to move to the nodes with greater degree (hub nodes), 
which enables the hub nodes to be efhciently used and 
enhance the system's overall capability; but when (j) is 
too large, the hub nodes will more probably get jammed, 
and the efficiency of the system will decrease. 

Finally, we study the fundamental diagram of net- 
work traffic with our model. Fundamental diagram (flux- 
density relation) is one of the most important criteria 
that evaluates the transit capacity for a traffic system. 
Obviously, if the nodes are not controlled with the queue 
length L, the network system will not have a maximum 
number of packets it can hold and the packet density can 
not be calculated, so that the fundamental diagram can 
not be reproduced. 

To simulate a conservative system, we count the num- 
ber of removed packets at each time step and add the 
same number of packets to the system at the next step. 
The flux is calculated as the number of successfully de- 
livered packets from node to node through links per step. 
In FiglH the fundamental diagrams for (j) = 0.0, 0.3, —0.5 
and —0.7 are shown. 

The curves of each diagram show four flow states: free 
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FIG. 3: (color online). Fundamental diagram for a, N — 1000 
network with mO = m = 5, a — 1, (3 — 0.2, and different 
(j). The data are averaged over 10 typical simulations on one 
realization of network. In each chart, the solid square line 
shows the flux variation when adding packets to the system 
(increase density), while the empty circle line shows the fiux 
variation when drawing out packet from the system (decrease 
density). The sudden transition density values are: 0.26 and 
0.23 (0 = 0.0), 0.40 and 0.34(0 = 0.3), 0.26 and 0.15((^ = 
—0.5), 0.15 and 0.13{4> = —0.7). For different realizations of 
network, the fundamental charts are similar, but with small 
difference in the transition values. The arrows in charts of 
4> = 0.3 and —0.5 are showing the hysteresis for guide of eyes. 



flow, saturate flow, bistable and jammed. For simplic- 
ity, we focus on the (j) = 0.3 chart with the maximum 
(Flux) = 1319 in the following description. As we can 
see, when the density is low (less than « 0.10), all pack- 
ets move freely and the flux increases linearly with packet 
density, which is attributed to a fact that in the free flow 
state, all nodes are operated below its maximum deliv- 
ering ability C. Then the flux's increment slows down 
and the flux gradually comes to saturation (0.10 ~ 0.34), 
where the flux is restricted mainly by the delivering abil- 
ity C of nodes. 

At the region of medium density, the model repro- 
duces an important character of traffic flow - "hysteresis" , 
which can be seen that two branches of the fundamen- 
tal diagram coexist between 0.34 and 0.40. The upper 
branch is calculated by adding packets to the system, 
while the lower branch is calculated by removing packets 
from a jammed state and allowing the system to relax 
after the intervention. In this way a hysteresis loop can 
be traced (arrows in Fig[3]), indicating that the system is 
bistable in a certain range of packet density. As we know 
so far, it is the first time that the hysteresis phenomenon 
is reported in the scale-free traffic system. 

One can also notice that when (j) = 0-3, the maxi- 
mum saturated (Flux) is higher than others, and the 
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FIG. 4: (color online), (a) Fundamental diagram for a. N = 
5000 network with mO = m — 5, a — 1, l3 = 0.2 and (f) = 
0.1. (b) The averaged number of jammed nodes (A'j„). The 
symbols for increasing/decreasing density are the same as in 
Fig|3] One can see that the two sudden change points 0.40 
and 0.14 in both charts are equal. The arrows are showing 
the hysteresis for guide of eyes. 



saturated flow region is much boarder than the cases of 
(j) = 0.0, —0.5 and —0.7. All these results show that the 
system can operate better when is slightly greater than 
zero, which is also in agreement with the simulation re- 
sult of Rc in Figdl 

In order to test the finite-size effect of our model, we 
simulate some systems with bigger size. The simulation 
shows similar phase transition and hysteresis in funda- 
mental diagram as shown in FigHI^a). 

The flux's sudden drop to a jammed state from a 
saturated flow indicates a first order phase transition, 
which can be explained by the sudden increment of full 
(jammed) nodes in the system (See FigHJ^b)). Accord- 
ing to the evolutionary rules, when a given node is full, 
packets in neighboring nodes can not get in the node. 
Thus, the packets may also accumulate on the neighbor- 
ing nodes and get jammed. This mechanism can trigger 
an avalanche across the system when the packet density 



is high. As shown in Fig^b), the number of full nodes 
increase suddenly at the same density where the flux drop 
to zero and almost no packet can reach its destination. 
As for the lower branch in the bistable state, starting 
from an initial jammed configuration, the system will 
have some jammed nodes that are difficult to dissipate. 
Clearly, these nodes will decrease the system efficiency 
by affecting the surrounding nodes until all nodes are 
not jammed, thus we get the lower branch of the loop. 

In conclusion, a new model for scale-free network traf- 
fic is proposed to consider the nodes' capacity and de- 
livering ability. In a systemic view of overall efficiency, 
the model reproduces several significant characteristics 
of network traffic, such as phase transition, and for the 
first time, the fundamental diagram for networked traffic 
system. Influenced by two factors of each node's capa- 
bility and navigation efficiency of packets, the optimal 
routing parameter is found to be slightly greater than 
zero to maximize the whole system's capacity. A special 
phenomenon - the "hysteresis" - is also reproduced in the 
typical fundamental diagram, indicating that the system 
is bistable in a certain range of packet density. It is the 
flrst time that the phenomenon is reported in networked 
traffic system. For different packet density, the system 
can self-organize to four different phases: free-fiow, sat- 
urated, bistable and jammed. 

Our study may be useful for evaluating the overall effi- 
ciency of networked traffic systems, and the results may 
also shed some light on alleviating the congestion of mod- 
ern technological networks. 
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